Radiation-condensation instability in a four-fluid dusty plasma 
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In this work linear stability analysis of a four-fluid optically thin plasma consisting of elec- 
trons, ions, neutral atoms, and charged dust particles is performed with respect to the radiation- 
condensation (RC) instability. The energy budget of the plasma involves the input from heating 
through photo-electron emission by dust particles under external ultraviolet radiation as well as 
radiative losses in inelastic electron-neutral, electron-ion, neutral-neutral collisions. It is shown that 
negatively charged particles stimulate the RC instability in the sense that the conditions for the 
instability to hold are wider than similar conditions in a single-fluid description. 
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I. INTRODUCTION 



The radiation-condensation (RC) instability in an op- 
tically thin plasma plays an important role both in space 
and laboratory plasmas (cf. Refs. [H, H, H, |3) • Dust impu- 
rities can alter dynamical and thermodynamical proper- 
ties of the plasma because of their high inertia and the 
ability to transform thermal energy into radiation. It 
was shown by Ref. within the framework of a two-fluid 
(dust and plasma) description that the presence of dust 
particles can change the conditions for the RC instability 
to grow. An essential assumption of Ref. was that dust 
particles contribute to the radiative cooling. In many 
situations, however, dust provides instead the energy in- 
put to plasma due to photo-emission of electrons from 
the dust surface when the system is exposed to ultravi- 
olet radiation. Such conditions can be met in interstel- 
lar plasmas (cf. Ref. 0). In this case the energy gained 
from non-thermal electrons escaping photo-ionized dust 
grains is transmitted to thermal electrons in the ambi- 
ent plasma, which then in turn heat ions and neutrals in 
elastic collisions. Simultaneously, the electrons lose their 
energy and pressure radiatively, and therefore are com- 
pressed under the external pressure. In this picture dust 
particles and electrons are the thermodynamically active 
components, while neutrals and ions represent the passive 
components providing conditions for radiative cooling of 
the electrons. Therefore, at such circumstances only a 
full four-fluid description is adequate. For the sake of 
simplicity we restrict ourself in this paper to considera- 
tion of a non-magnetized plasma. 

The paper is organized in the following way: in 
Sec. II. A we present the full set of four-fluid equations 
describing the evolution of the plasma under considera- 
tion, followed by the description of the energy exchange 
between the fluids in Sec. II. B. In Sec. III. A we present 
the linearized dynamical equations, and in Sec. III.B we 
derive the instability criterion for the condensation mode. 
The analysis of the instability is given in Sec. IV. We close 



our paper with a short summary in Sec. V. 



II. DYNAMICS AND THERMODYNAMICS 
A. Dynamical equations 

We start from a four-fluid unmagnetized system (based 
on Refs. and d for a magnetized plasma) for electrons 
(a=e), ions (a=i), neutral atoms {a=a) and massive 
charged dust particles {a=d): 



dtUa + V • {UaUa) = 0, 



(1) 



Vpa V$, (2 



dtPa = -V • (paUa) + (7 " 1)(-PqV • + Ha), (3) 



(4) 



where Ua, Za, and Pa are particle number density, 
electrical charge number, velocity and gas pressure of the 
ath species, respectively. As shown in Ref. for a wide 
range of wavelengths relevant to the RC instability, sep- 
aration of the dust charge can be neglected, such that 
the charged components are coupled through the quasi- 
neutrality condition 



(5) 



We assume the ions to be singly charged, so that 
Zi = +1 and = — 1. As the dust is usually neg- 
atively charged in the systems under consideration, we 
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write Zd = —Zd, so that a positive value of Zd corre- 
sponds to negatively charged dust and vice versa. The 
quasi-neutrality condition, thus, reads 

rii - ZdUd - ng = 0. (6) 

However, the validity of this assumption is limited by 
the absolute value of dust charge: the characteristic time 
scale for a radiatively cooling plasma is of the order of 
^ ksT/ Kn ^ lO^^n^^s in temperature range T ^ 10'* K. 
Here, A ~ 10"^'^ erg cm^ s""'^ is the radiative cooling 
function at T ~ 10'* K for the interstellar environment, 
and n is the total number density of all species together. 
Dust separation becomes important when this time is 
shorter than the inverse dust plasma frequency uj^^ — 

^m.d/'^T^Z'^e^nd ^ lO^n^^/^ s, for iid ^ 10~*^n typical 
for interstellar plasma, which gives Zd ^ lQ~^y/n. 

For simplicity we neglect here the collisional source 
terms in the momentum equations. In the case when 
all the four fluids in the unperturbed state are at rest, 
Uq = 0, we consider here, this assumption can only result 
in an overestimate of the growth rate, but cannot change 
the instability criterion. 



B. Energy exchange 

The source terms in the energy equations, TLa = 
Ha — La, describe the generalized energy gain rate due 
to radiation and collisional processes. In the interstellar 
medium (ISM) the RC instability develops mostly in the 
warm neutral medium with T ^ 10* K, n ~ 0.1 cm~'^ and 
a fractional ionization of cc 0.1. Small amount of trace 
heavy elements account for ~ 0.01 of the gas mass, and 
a similar amount is confined in dust grains. The sizes of 
grains range from ~ lOA to 0.4/Lim (cf. Ref.Q), the charge 
of dust grains in the warm neutral medium is positive 
and varies from ~ 0.1 to ~ 10^, e.g. Ref. Hoi, 
depending on size; for the sake of simplicity we will de- 
scribe in what follows the dust component as an ensem- 
ble of particles of equal masses, md — 10~^'*g, and equal 
charges, Zd, assuming the average abundance of dust par- 
ticles, Xd — nd/n ~ 10~^^. The ion component comprises 
mostly of protons and singly ionized trace elements such 
as CII, Fell, and Sill (ionized by the interstellar ultra- 
violet field), the neutrals contain hydrogen atoms, HI, 
and neutral oxygen, 01. Helium is neutral in the warm 
neutral medium, and its contribution to radiative pro- 
cesses is negligible at low temperatures T < 2 x 10*^ K, 
therefore it can be included in the system by multiplying 
the mass of the neutral particles by factor of 1.4. We 
have used the two sets of abundances corresponding to 
the standard solar abundances Xi = rii/n for species i: 
XFe = 3.2 X 10-^ xs; = 3.2 X 10-^ xq = 4.4 x 10"'*, and 
xc — 3.57 X 10^*^, and a "depleted" set of abundances 

with XFe = 6.0 X 10"'^, XS^ = 2.0 X 10-^ Xq = 5.0 X IQ-*, 

and Xc = 6.0 x 10~^, describing the mean composition 



of the interstellar plasma with some elements frozen on 
dust grains as suggested first by Ref. JJj. The results for 
the two sets are similar; later only calculations for the 
"depleted" set are shown. 

We assume in this paper that the system gains the en- 
ergy from photo-electrons produced by external ultravi- 
olet radiation ionizing dust particles. Non-thermal pho- 
toelectrons emitted by dust share their energy with ther- 
mal electrons of the plasma with the rate Fn^, where 
F is determined by the external UV radiation flux. Go, 
and optical properties of the dust. In inelastic collisions 
with ions and neutrals the electrons lose their energy to 
excite internal degrees of freedom, which then decay ra- 
diatively. In the simplest case of an optically thin low- 
density plasma these energy losses by the electrons can, 
thus, be written as the sum of two processes: L^{Te)nine 
and L'l{Ti,)nane, where we implicitly assumed that the 
thermal velocity of the electrons VT.e 3> VT,i, VT,a, so 
that in LI and only the dependence on is essential. 
In addition, the electrons can share their thermal energy 
in elastic collisions with the ions and neutrals. The coef- 
ficients are and q'^, respectively. In total it gives 

He = Tud - L^(Te)nine - Ll{Te)nane 

- ql{Te){Te - Ti^riiUe 

- ql{Te)(Te -Ta^naUe, (7) 

where for Tud we have taken the value typical for the 
interstellar plasma illuminated by the galactic ultraviolet 
(UV) radiation field according to Ref. \1% The cooling 
functions Lf and ijj are taken from Refs. Il2l and [H, the 
coefficients for the collisional energy exchange between 
the electrons and ions qf = qf{Te) and the electrons and 
neutrals q% = qa{Te) are given in Refs.[l3 andlTsl. 

Energy losses of the ions in inelastic collisions are small 
in comparison with energy exchange rates in elastic col- 
lisions, and these losses are normally neglected. In our 
calculations we included, however, these losses, assum- 
ing, according to Ref. 16, that the cross section for direct 
excitation by proton impact is broadly similar to the one 
by electron impact if the proton has the same velocity 
as the electron, with the only difference that the max- 
imum cross section by proton impacts is a factor of 10 
larger than that by electrons (cf. Refs.[l3,[ll[li). With 
this assumption the derived form of the radiative cooling 
functions for the ions read as L\ ~ a{me/miY^^L': for 
the ion-ion inelastic collisions and ~ a(me/TOi)^/^L^ 
for the collisions between the ions and neutrals, where 
we varied the factor a from 1 to 10. Therefore, we adopt 
for the ions the generalized rate in the form 

Hi = -L\{Ti)ni^ - L\{Ti)nani - qd{Ti)Tinind 
-qt{Te) {Ti - niUe 

-q'a{Ta,Ti){T,-Ta)n,na, (8) 
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wher e q!, de pen ds on both Ti and Ta and is given in 
Refs. IlirandflSl: describes coohng of the ions in elas- 
tic coUisions with translationally cold dust particles - by 
the order of magnitude 5^ ^ 2kB{'mp/md)TT(Td^VTi, where 
TT^p/'m-d is the mass ratio of the ions (mainly protons) to 
dust particles, (Jd the geometric cross-section of a dust 
particle, the thermal velocity of the ions, and is small 
compared to other terms in ([5]); the second term on the 
r.h.s. describes radiative losses of the ions in collisions 
with neutrals; the contribution of ion- ion inelastic colli- 
sions in a relatively weakly ionized plasma we deal with 
{xe ^ 0.1) is small. 

Similarly we obtain for the neutrals 



• Q • , , 1 ",3,0 
— tUJiya — U^t>a = —l^KUa + / ^ — — — 



„ Ta,0 dnp Ua^O 



E 



1 dhg Tp^o 

nafi dTp Tafi 



47re \ - ~ 



(14) 



where = (7 - l)Tia/kB, = Sna/uafl, da = 
STa/Ta,o, ks is the Boltzmann constant. Combining 
(fTTj) . p^ . and (HSl) one obtains 



Ha -LfiTa)n,na~L''^iTa)nl~q^{Ta)Tanana 

-ql{T,){Ta-Te)nane, (9) 

where the first term in the r.h.s. describes radiative en- 
ergy losses of the neutral atoms (mostly hydrogen) in 
collisions with the ions (mostly the ions of heavy ele- 
ments, like CII, Fell), while the second term corresponds 
to energy losses in collisions of mostly hydrogen atoms 
with such heavy elements as 01. The third term in ^ 
describes cooling of the neutrals in elastic collisions with 
presumably cold dust particles - loosely ~ 9d ^^"^ both 
contributions are small in comparison with other terms 
and are, thus neglected in the results shown below. 

In the equilibrium state of the component the following 
set of equations are fulfilled 



Tie = 0, TLi = 0, Tia = 0, 



(10) 



where the energy gain for the ions and neutrals stems 
from a small difference in temperatures Tg > T^, Ta, such 
that Te — Ti, Te — Ta <^ Te and approximate equalities 
Te Ti Ta hold in the whole temperature range. The 
solution of (fTU]) is depicted as dotted lines on Figs 1-3. 



III. STABILITY ANALYSIS 

A. Linear perturbations 

We consider planar perturbations, Sg, of a quantity g 
in the form Sg cx exp[i{kx — i^t)], i.e. we assume the per- 
turbations to propagate only in x-direction for simplicity. 
In this case the equations ([1]) to ([4]) read 



UJUc 



iOVa — kUa = 0, 



kc^Va + kCa-Oa + 



i(7 - l)ujVa 



iujriap + \al3 + ^ A^qk?7k/3 ) = 0, 



(15) 

where a and k, run over the lighter species i, e, and a, 
whereas index (3 runs over all four species, i.e. i, e, a, and 
d. The matrix rjap is given in Appendix A. This equation 
must be complemented by a condition that follows from 
inserting Eq. dH]) into Eq. for the (cold) dust: 



'e 2 



rrid 



Zd—ujp^Ui + {uj d )vd^ 0, (16) 



where ojpa is the plasma frequency of the ath component. 
Further we define the abbreviations 



A 



1 dhg npfl 

Tafl dnp Uafi ' 



1 dhaT 



/9,0 



"a,0 Tafi 



(17) 



Explicit expressions for Aq^ and flap are given in Ap- 
pendix B. We concentrate in this paper on the case when 
the external (unperturbed) electric and magnetic fields 
are zero, and, thus, in linear approximation the Lorentz 
force can be neglected. 



B. The instability criterion 



From system ([T5|) and Eq. p6)) we obtain the following 
dispersion, where we use the abbreviation N — —ilu for 
the growth rate: 



-rjeiN+Aei 

Aai 
nij y 2 



-?7,eiV-fA,e A 
L'eA^+Aee A, 

^^Zdule 



rf,dN+K.M 

r]edN+Aed 
aa ^ad 



with A, 



Q/3 



Aq/3 + Z^KA^aK^ft/S' -Da = (7 - 1) 



(18) 



In the limit <C iLJpd ^^'^ charge separation becomes 



(11) unimportant and converges to 



(12) 



D^N + A„; 


A,e 


Ka 


Aid 


Aee 


DeN + Aee 


Aea 


Aed 


Ki 


Aae 


DaN + Aaa 


Aad 


-n^ZdUd 


ne/ZdUd 





1 



= 0, (19) 
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where K^p = Xaf3 + Hafidfj, dp = a;^/fc^c| — 1. We as- 
sumed here exphcitly that the heating and coohng rates 
of the plasma components e, z, a do not depend on the 
dust kinetic temperature, T^, which impUes pL^d — IMd = 
fJ'ad = 0. Equation (fT9)) is the dispersion relation for a 
radiating plasma, when quasi-neutrality can be assumed. 

In the low- frequency limit |A'^| <C kce, kci, kca the 
coefficients Da converge to Da = 7, and therefore the 
sufficient condition for the instability is that the constant 
term in the polynomial in N is negative. Thus, the 
instability condition is 



ZdTld 



rii 
Zdfid 



XadKiKa 




+ ^edA'uKa 


XedJ^'aiKa 


\ A ' A ' 


- Xidh-'^iA'a^ 


\ A ' A ' 

^ad^^ee'-^ia 


+ ^idA^gAgjj 


+ Xed^'ie^'aa 


XedKa^'ae 


\ A' \' 


- Xid^-'ee^'aa^ 


A' A' A' 

ee 11 aa 


+ A' A'-A' 

' ae ei la 


+ A' A' A' ■ 

' le ea at 


A' A' A' 

ae ea 11 


- A' A' A' 

le ei aa 


- A' A' A' 1 

ai la ee j 



<o, 



(20) 



here A' 



Q/3 



Xa(3 — fJ-a/n- This Condition is an extension of 
the standard (one-fluid) condition for the condensation 
mode to grow as given by Ref . H: 



dH To m 



IV. RESULTS 



(21) 



In a one-fluid description under the conditions consid- 
ered here: H oa nd oc n and L = A{T)n? with fixed 
(temperature independent) ionization states, the insta- 
bility criterion reads as dlnA/dhiT < 1 and restricts a 
density independent temperature range where the plasma 
is unstable: for the diffuse warm neutral phase of the 
ISM the instability can occur in the temperature range 
100 ;S T ^ 7000 K. The situation changes dramatically, 
when a possible separation in the dynamics of the four 
plasma components is taken explicitly into account. For 
the sake of simplicity we consider first constant dust heat- 
ing rate r=const in the whole range of temperature and 
density. This suggests that the dust charge and corre- 
spondingly the fraction of the UV radiation transformed 
into gas heating are kept constant in the whole range of 
plasma parameters. It should be stressed though that at 
realistic conditions in the interstellar plasma dust charge, 
and as a result the heating rate, depend on temperature 
and density, cf. Ref. 0: Zd = Zd{T,ne), T = r(T,ne). 
However, we believe it worthwhile to restrict ourselves 
with a simplified scheme in order to understand bet- 
ter how dust particles affect the dynamics of the RC 
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FIG. 1: The instability criterion in the temperature-density 
plane in a double logarithmic representation. Grey areas de- 
pict unstable regions, white corresponds to stable perturba- 
tions; the dotted line is the solution of the equilibrium con- 
dition r = A(r)n written in the one-fluid approximation and 
describes the equilibrium equation of state T = T(n). Panel 
a) shows the separation of the temperature-density plane into 
stable and unstable domains for the standard one-fluid case), 
panel b) depicts the condition for Zd = Q (b). The middle 
row shows the shows the separation for positive dust with 
Zd = —1 (c) and Zd = —10 (d). The bottom row shows the 
cases Zd ~ —100 (e) and for negative {Zd = 0.1) dust (f). 
Note that all cases for negative dust look almost identical in 
the plane shown here. 



instability in a four-fluid description. In the models 
shown below we assume also a constant fractional ion- 
ization, Xe =const. This assumption is justified by the 
fact that the instability growth time is expected to be of 
the order of the radiative cooling time, which is always 
shorter than the relaxation time for recombination (see 
Ref. [20h ; for conditions, e.g., in the warm phase of the 
interstellar medium with T 10'* K the radiative cooling 
s"""^, while the recombination time is 



10 



12^ 



ill, 

-1 



Figure 1 shows the domains of instability in the 
temperature-density plane for several values of the dust 
charge. Panel a) corresponds to the standard one-fluid 
case for a constant fractional ionization = ne/n=0.1: 
the whole plane is separated onto two regions inde- 
pendent of the density, where the plasma is stable 
(white) or unstable (gray), corresponding to the criterion 
d\nA/d\nT < 1. Panel b) shows the instability domain 
in the four-fluid description with neutral dust particles 
Zd — 0: the domain is clearly density- independent and 
slightly wider than on the panel for the one-fluid case 
- this corresponds to the condition dlnA/dlnT < 2, 
which holds either when radiative losses are unbalanced 
by heating, or when heating is independent of the den- 
sity, i.e. H =const. As we neglect in this study collisional 
friction between the components, the dust component is 
decoupled from the plasma when Zd — 0, which is equiv- 
alent to the latter case, i.e. SH — 0. 

The panels c), d), and e) correspond to positively 
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FIG. 2: The instability criterion on the temperature-density 
plane for highly charged negative dust grains; Zd ~ 10^. 

charged dust grains, Zd < with increasing absolute 
values (given in the Figure caption) of the dust charge. 
It is clearly seen that positively charged dust particles 
destabilize the system, widening the instability domain 
to higher temperatures where formally in one-fluid ap- 
proximation dlnA/dlnT > 1. 

Negative dust particles, shown in panel f), strongly 
stabilize the system even in the limit of small Zd - dust 
charges in the range Zd = +[0.01 — 100] give identically 
stability in the whole domain of the temperature-density 
plane shown in Fig. 1. Stabilization of the thermal in- 
stability by negatively charged dust can be understood 
as follows: due to the quasineutrality condition negative 
dust particles repel the electrons out of the compressed 
plasma, and thus the radiative cooling rate provided basi- 
cally by the electrons decreases. Only when dust charge is 
high in absolute value, Zd = 10'^, a narrow region at tem- 
peratures 100 K ^ K, where the instability con- 
dition is fulfilled, does appear in the high density range 
as seen in Figure 2. 

This instability at high densities in turn seems to stem 
from a tighter collisional coupling between the electrons 
and ions. Figure 3 shows the dependence of the insta- 
bility domains in the temperature-density plane for an 
artificially enhanced rate of the elastic energy exchange 
between electrons and ions; the left panel corresponds to 
positively charged grains, the right panel to negatively 
charged grains; in both cases the dust charge is fixed at 
the level jZ^I = 10'^. It is readily seen that weakening 
of the elastic energy exchange in collisions between the 
electrons and ions results in a widening of the instability 
domain for positively charged grains (left upper panel), 
and in a shrinking of the instability domain for nega- 
tively charged grains (right upper panel). More frequent 
electron-ion collisions lead to a narrower instability do- 
main for positive dust, while in the case of negative dust 
the instability domain widens with the collisional rate. In 
the limit of very high collisional frequencies (increased by 
a factor of 100 with respect to the normal value) the in- 
stability domains for positive and negative grains become 
almost identical. This means that the instability domain 
at high densities for negatively charged grains shown in 
Figure 2 arises due to elastic interactions between the 
electrons and ions. The dependence of the instability do- 
main on the density emerges from the fact that the energy 
exchange in elastic collisions between the species does not 
appear in the equilibrium energy equation because of the 
isothermality of the equilibrium Te = Ti = Ta, while in 
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FIG. 3: The instability criterion in the temperature-density 
plane for enhanced rate of elastic energy exchange between 
electrons and ions; the left and right panels show the in- 
stability domain for positively {Zd = —10^) and negatively 
{Zd = 10^) charged grains, respectively. The collisional fre- 
quency is multiplied by a factor of 0.01 (a,b), 10 (c,d), 20 
(e,f), 50 (g,h), 100 (i,j). 

the instability criterion ([20)1 the derivatives of the elastic 
collisional rates over temperatures, proportional to the 
density square, do contribute; at high densities this con- 
tribution becomes dominant. 

For a constant dust charge Zd variations of the UV ra- 
diation flux. Go , result in a simple shift of the instability 
domain and the equilibrium temperature curve along the 
n-axis - this is a consequence of the fact that at such an 
assumption the heating rate linearly depends on density 
n, while the cooling rate is proportional to the square 
of density n^: the higher is Go the larger is the density 
needed for radiation losses to balance the heating at given 
temperature. This behavior may change to some extent 
when the dependence of dust charge on Go is accounted. 



V. SUMMARY 

The radiation-condensation (thermal) instability in a 
dusty plasma heated through photo-ionization of dust 
particles is described in a four-fluid approximation. It is 
shown that in this approximation with admitted separate 
motions of the species, the instability criterion changes 
dramatically compared to the standard one-fluid case. 
In particular, positively charged dust particles are found 
to destabilize the system, while negatively charged ones 
strongly stabilize perturbations: only in the limit Zd^ 1 
negative dust leaves a relatively narrow domain of insta- 
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bility in the high density end, where elastic energy ex- 
change between the electrons and ions dominates. 

In the interstellar medium the dust charge in low- 
density regions (diffuse neutral HI phase) is positive due 
to photo- ionization by stellar ultraviolet light, cf. Ref . ITol. 
In the light of our results this means that low-density in- 
tercloud gas is highly unstable against the formation of 
condensations in the whole temperature range. A signif- 
icant fraction of dust is also positive in denser regions, 
such as diffuse HI clouds and even in molecular clouds 
(see Ref. [To| ). Therefore, contrary to a common under- 
standing that dense phases of the interstellar gas are sta- 
ble against the formation of condensations, they turn out 
to be unstable when separated motions of the species are 
accounted. One may speak thus about possible fragmen- 
tation of interstellar clouds through the RC instability 
on smaller structures. 



T 

^ e. 

■^ed — — 7^ 



Ana = ^ 





+ 2Llna 



(B2) 



Ta 



Kd = -qd'^d 



(B3) 



APPENDIX A: r^^p 

The matrix rjap is determined from the momentum 
equations for the components and connects the dimcn- 
sionless temperature and density 



= ^ riapvp 

13 

The elements of the matrix are 



f IJii 1]ie f}ia f]id 

f]ei f]ee 1]ea '^ed 

^ai Vae ^]aa ^ad 

\ Vdi Vde Vda Vdd 

,22 



(Al) 



mi k^c? 








i 






(A2) 



APPENDIX B: A^a AND 

For simplicity we assume further that F does not de- 
pend on Tie and T^. In addition, we assume that in equi- 
librium |Te - r,|, iTe-Tal, \T,-Ta\<^T^. The collision 
coefficients with the dust, 9^ and give here, were ne- 
glected in the results shown. The coefficients are: 



\i — 

\e = 

Aia — 

Aid = 



2L\n, + 


T * 
^ i 

^Qd^d 



Qd^d 



^J'id 



dL\ dLi dql, 

rii Tin nrf 

dT^ dT, dT^ 

-itne - qlna - qdUd, 



= 



(B4) 



dLf 



dL't 



dTf, dTf. 
-lint - qlna, 





(B5) 



fJ-ad 



dLf dLl dq2 

-Ui -Ua —Ud 

dia dla dla 



0. 



(B6) 



(Bl) 



We accounted here that direct injection of energy from 
dust particles to the ions and neutrals is zero. 
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